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Abstract. The residence time of dust grains inside tur- 
bulent molecular clouds is calculated . Earlier estimates of 
this time by Boland & de Jong (1982) have been criticized 
by Prasad et al. (1987) for not taking into account prop- 
erly the turbulent character of the velocity field in molecu- 
lar clouds which implies that the trajectory of a dust grain 
should be treated as a random walk. Taking some minimal 
assumptions regarding the turbulent velocity field a time 
scale is derived that depends on the Reynolds number of 
the flow. The near proportionality with Reynolds number 
of this time scale results in a much longer time scale over 
which dust grains will remain inside a molecular cloud. 
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1. Introduction 

Chemical surface reactions on dust grains form an es- 
sential pathway in the formation of complex molecules, 
and even some quite simple ones, in interstellar molecular 
clouds. The time that a grain resides inside clouds with a 
high optical depth determines its life time and therefore 
the time over which it can play a role in the chemical pro- 
cessing of its host molecular cloud. This residence time 
scale has been estimated previously by Boland & de Jong 
(1982) to be L/a v from the observables L which is the size 
of the cloud and o~ v which is its velocity dispersion. Prasad 
et al. (1987) argued that this essentially implies a linear 
path through the cloud for any dust grain. This does not 
properly account for the fact that the velocity field in the 
cloud is turbulent and that the grains, dragged along by 
the gas motions, perform a random walk. Using a con- 
stant step size random walk Prasad et al. (1987) derive 
a time scale L 2 /a v l s where l s is the (typical) step size 
in the random walk, which they estimate to be the scale 
of the largest eddies in the cloud. However using a single 
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scale eddy size to represent the large range of eddy sizes 
present in highly turbulent flow is also inappropriate. Tak- 
ing dust grains to be test particles of the turbulent flow, 
it is possible to derive a circulation time scale in terms 
of the Reynolds number of the flow. From the study of 
Weidenschilling and Ruzmaikina (1994) it is clear that in 
particular sonic turbulence strongly affects grain growth 
and destruction since it affects the rate of collisions and 
the energy of such collisions with other grains and with 
gas phase material. The point of focus of this paper is 
not the growth and destruction of grains but rather the 
stochastic migration of an ensemble of grains throughout 
a dense cloud. It is demonstrated that for quite general 
properties of the turbulence the random walk through a 
cloud indeed considerably increases the time grains stay 
in a region of high optical depth, in accordance with the 
argument of Prasad et al. (1987) although the time scale 
derived here is different from theirs. 

2. Power law turbulence 

The coupling between gas and dust grains in stellar winds 
and in the interstellar medium can be shown to be quite 
strong, because of the strong drag that dust grains ex- 
perience when drifting with respect to the gas. The dust 
therefore is treated here as completely momentum coupled 
to the gas : effectively it is a test particle for the turbulent 
flow. 

The length scale of dissipation (internal scale) of the 
turbulence Id is set by the viscosity in the gas of the molec- 
ular cloud or the magnetic diffusivity (the inverse of the 
conductivity) if the cloud contains a (tangled) magnetic 
field. Unless this length scale of dissipation of the turbu- 
lent velocity field is nearly equal to the size of the molec- 
ular cloud the energy density E in the velocity field per 
length scale interval dl, or equivalently per wave vector 
interval dfc, will have an 'inertial range' over which it be- 
haves as a power law as a function of length scale. Outside 
of this range the energy density or power in the turbulent 
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velocity field decreases rapidly, cither due to dissipation or 
because the cloud is not large enough to contain the length 
scales larger than its own size. For stationary hydrody- 
namical turbulence in the absence of magnetic fields, and 
in the inertial range, the energy density has a slope of this 
power law which conforms to the well known Kolmogorov- 
Obukhov scaling : E(k)dk oc fc~ 5//3 dfc (cf. Landau & Lif- 
shitz, f 987). As was pointed out by Scalo (1987) molecular 
clouds are highly magnetized and subject to self-gravity. 
Also energy is injected into the turbulent velocity field 
from a variety of sources on a variety of spatial scales, 
for instance shear due to galactic rotation, so the slope 
of the energy spectrum can differ from the value of —5/3. 
However even for magnetohydrodynamic (MHD) turbu- 
lence the velocity field behaves as a power law over a 
large range of length scales (cf. Stanisic, 1985). Gravity 
does not impose a characteristic length scale nor do ef- 
fects of compressibility. The energy injection mechanisms 
can be subject to one or more characteristic length scales 
which could introduce humps or roll-overs in the turbu- 
lent spectrum, where locally the energy spectrum is not 
a power law. However the requirement of stationarity of 
the turbulent spectrum imposes strong constraints on the 
transfer of energy between scales. Only if the turbulent 
cascading of energy is the least effective (slowest) of the 
energy transfer mechanisms operating, which is unlikely, 
will the spectrum be severely affected by the other effects, 
because then energy can build up at the injection length 
scale before the energy cascading from larger scales and to 
smaller scales can compensate for this. Summarizing, the 
energy density probably does not follow exactly a power 
law with a slope equal to —5/3 over the entire range of 
scales possible, but the departures from this behaviour are 
not likely to be large. 

Obscrvationally it is relatively straightforward to de- 
termine a power law spectrum for the mass M contained 
in clumps within molecular clouds, so the starting point is 
this mass spectrum, where f(M) is the fraction of cloud 
clumps with a mass between M and M + dM : 



/(M)dM oc M- a dM 



(1) 



Recent measurements of a by Williams ct al. (f 995) yield 
a value for a = 1.27 ± 0.09 for the Rosette Nebula. Val- 
ues of up to 1.7 for this and other cloud complexes have 
also been reported (cf. Blitz, f99f). It appears that a is 
roughly constant over a range of scales from 100 pc down 
to 10~ 4 pc in at least one cloud (cf. Rouan et al. 1997). 

2.1. In the absence of density-velocity correlation 

The mass is related to the gas density p and I by : 

M oc pi 3 (2) 

In principle the gas density could be a function of the 
length scale I, which is discussed in Sect. 2.2. Assuming a 



p independent of I the mass spectrum can be converted to 
a spectrum in terms of sizes between I and I + dl : 



f(l)dl 



c\i 2 - 3a di 



(3) 



Outside of this range the amount of material is assumed to 
be negligible and the value of the power-law index a is as- 
sumed constant over the entire inertial range. Even though 
the observed clumps should probably not be individually 
identified with turbulent eddies it is likely that the turbu- 
lent motions follow a similar distribution of energy density 
versus size, since the clumps are a consequence of random 
convergence of the turbulent flow. In this paper it is as- 
sumed that the distributions arc identical. The constant 
Cl is determined by demanding that the integral of /(/) 
over I between the dissipation length scale Id and the size 
of the molecular cloud L be equal to unity. 



C L l 2 ~ 3a dl = 1 



implies : 



(4) 
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9a-9 

Re 4 "I 



(5) 
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a = 1 
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where the Reynolds number R e is related to L and Id by : 

4/3 



" = (0 



(6) 



If a cloud is optically thin in the spectral line in which 
the velocity dispersion is measured the observed velocity 
dispersion of a cloud must be the expectation value of the 
square of the velocity. If the cloud is optically thick the line 
only measures the turbulence in the surface layers. Since 
the turbulent velocity field is probably not homogeneous a 
velocity dispersion determined with an optically thick line 
is then not representative of the interior of the cloud. In 
the following it will be assumed that a measured velocity 
dispersion is representative however, so that either the 
method with which it is measured is sensitive to velocity 
fields throughout the clouds or the turbulence is assumed 
to be homogeneous and isotropic. In the usual turbulent 
cascade (cf. Landau & Lifshitz, 1987) of fully developed 
free turbulence the velocity vi of a turbulent eddy with 



size f is 



vi 



VL 



L 



1/3 



(7) 



where vl is the overturning velocity of the largest eddies 
supported by the molecular cloud which have a size of the 
order of the cloud size L. As argued at the start of this 
section a number of effects could change the power I /3 in 
relation (7), but the departures are not likely to be very 
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large. The velocity dispersion a 2 is the expectation value 
of the square of the turbulent velocity : 



°l = J C L vfl 



dl 



C L v 2 L L 3 - 3a 
9a - 9 



11 -9a V 



l-R'e 



9a - 11 



11 (8) 
9 
11 



In the approximate equality it has been assumed that the 
Reynolds number R e ^> 1. The exponent 7 satisfies : 



7 = 



-1/2 

9-9a 
4 





a>f 
Ka<f 
a < 1 



(9) 



Note that the intrinsic thermal width of the line is as- 
sumed to be much smaller than the width of the velocity 
distribution, which implies that the turbulent velocities 
are supersonic. If this is not the case a proper convolu- 
tion should be done of the line profile and the velocity 
distribution. 

As a grain is swept along in the turbulent eddies it 
performs a random walk through the cloud. If all the steps 
were of equal length I then it would take Ni — (L/l) 2 steps 
to traverse a distance L, each step taking a time Ti = l/vi. 
The step size obeys the distribution (3) and so the total 
residence time for a dust grain must be the expectation 
value of NixTi : 



T tot = J ClNiTiI 2 -^ 

L 4- 3a 3 

= W 



vl |5 — 9a I 
(9a - 9) 3 



(1-iC*- 5 ) 

1/2 



a^ 



9 (10) 



(9a- ll)(9a-5) z 



tdS — ft L r>S 
n e = C T ti e 
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where the exponent S satisfies : 



3/4 
(17-9a)/8 
(9a-5)/4 




a > 
1 < a < 



§ < a < 1 
a < I 



11 
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(11) 



The slope a of the mass spectrum of the clouds appears 



in the exponent S only for 



< a < 



11 
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For a = 1.27 



the exponent S = 3/4 and Ct ~ 0.9. For steeper (a > 
1.27) spectra Ct approaches unity asymptotically. It is 
only for mass spectra that are flat or rising (a < 0) that 
the estimate of Boland & de Jong (1982) is appropriate, 
because only then is 8 = 0. 



Of interest is the dependence of the residence time T tot 
on the Reynolds number R e for a > 5/9. It seems perhaps 
counter-intuitive that for more vigorous turbulence (larger 
R e ) dust grains take longer to traverse a cloud. However 
one should realize that if the turbulence is more vigorous 
the velocity dispersion <j v also increases. For a given (mea- 
sured) velocity dispersion an increased Reynolds number 
implies that Id is smaller so that there is a larger contribu- 
tion to the integrals (8) and (10) from small scales. Thus a 
grain is likely to make many more small steps in a random 
walk than large ones. 

2. 2. Fully coupled density-velocity fluctuations 

In Sect. 2.1 the density p is assumed to be constant, so that 
it can be ignored in the integral for the velocity disper- 
sion. In view of the Larson relations (Larson, 1981) which 
show that the density of clumps follows a power law as a 
function of clump size, it is instructive to consider also the 
effect of a gas density which fluctuates in a manner that is 
spatially correlated with the turbulent velocity and which 
has a power law spectrum of I : p cx The measured 
value of (3 = 1.3 ±0.2 (cf. Cernicharo, 1991 and references 
therein). The combined distribution function is assumed 
to be the product of the individual distribution functions 
for the density and the velocity, which implies that (3) 
becomes : 



f(l)dl cx l 2 - 3a+ ^ a -Vdl 



(12) 



In this paper the clouds are not assumed to be in virial 
equilibrium, thus no a-priori relation between pi and vi is 
postulated, and j3 is a free parameter. For a cloud that 
is in virial equilibrium on every scale one expects that 
the kinetic and potential energies of turbulent eddies are 
identical up to a constant factor. This implies that vf cx 
Mi/l oc pi 2 oc / 2_/3 . Combining this with (7) shows that 
then one would expect (3 = 4/3, which is equal to the 
value actually observed to within lcr. 

In the analysis presented here a v depends on p since 
the density plays a role in the determining the profile of 
the spectral line from which the velocity dispersion is mea- 
sured. If one takes a density weighted average instead of 
(8) then : 



L, 

7)! ClPiV ' 



2;2- 3a+/ 3(«-l) d ; 



(P) 

Id 

I i ? (9«-9-3/3(a-2))/4 ( 13 ) 

L (p) L ll-9a + 3/3(a-2) 



a^ 



6/3-11 



3/3-9 

where (p) is the expectation value for the density. 
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Fig. 1. Ttot is a power law of R e for R e S> 1. The sign of 
the power law index is shown as a function of the power-law 
exponents a and (3. It is positive in areas indicated by a +, 
negative where there is a — and zero on the solid curves and in 
areas indicated by 0. On the dashed curves T to t oc (fn_R e )~ 1//2 . 
On the dash-dotted curve T to t oc In R e . The error bars show the 
range of measured values for a and (3 for clumps in molecular 
clouds 

Note that there is a problem with the velocity disper- 
sion that results from this analysis. Observations appear 
to indicate that <r v oc L 3 (cf. Scalo, 1987), whereas with 
the observed values of a and (3 one obtains a v oc v\R e 
which is independent of L because of relations (6) and (7). 
One way to explain this is that when observing at increas- 
ing resolution (smaller L) one concentrates on 'clumps', 
which are the regions with higher density and lower veloc- 
ity dispersion. This 'selection effect' might introduce the 
observed trend, bringing it closer to relation (7). 

The relation (10) becomes now : 



Ttot = J c L Nmi 2 - 3a di 
Id 

T n 1 R (9a-5-3/3(a-l))/4 ( 15 ) 

a v v L L 5 - 9a + 3/3(a - 1) 



a? 



3/3-5 
3/3-9 



It is not particularly illuminating to show the equations 
for all possible values of the exponents a and (3, since 
this separates out into 73 separate cases. Instead in Fig. 
1 is shown the behaviour of the exponent of the Reynolds 
number R e as a function of a and (3. Indicated by +, 0, 
and — are regions where this exponent is positive, zero, 
or negative respectively. For 2/3 < /3 < 3 and a > 1 
and for j3 < 2/3 and a > 2 the dependence of T tot on 
the Reynolds number is always Ttot oc for R e 1. 
From Fig. 1 it is clear that only for quite steep distribu- 
tions for the density (large (3) the exponent is negative and 
Ttot decreases as R e increases. In such cases the velocity 



dispersion is weighted strongly towards the high density, 
low velocity clumps so that for more vigorous turbulence 
the estimate T to t = L/a v is increasingly an overestimate 
of the residence time. For the measured a = 1.27, T to t 
is a decreasing function of R e only for (3 > 6.7. For the 
observed range of values a and (3 the relation is 



L L 2 

Ttot — Ct — R^ 4 = Ct-, 

a v Ldcr v 



(16) 



where the value of Ct ~ 0.35. This appears to be very 
similar to the expression of Prasad et al. (1987) who have 
l s replacing l d , where Z s is the typical eddy size. The dis- 
sipation length scale Id is much smaller than the typical 
eddy size I4 <C l s however, so that the time scale derived 
here is longer than that derived by Prasad et al. (1987). 

2.3. The influence of fractal structures 

It is important to note that the residence time in the cloud 
complex derived here may be an over-estimate of the res- 
idence time of the grains. There is growing evidence that 
the gas density distribution in molecular cloud complexes 
is highly inhomogeneous or 'fractal' in structure (cf. Scalo, 
1987, 1990; Falgarone et al., 1991). Even though the grain 
can take a long time to traverse a cloud complex it will 
encounter regions within that complex that may well be 
sufficiently warm or exposed to ultraviolet (UV) radia- 
tion to be quite destructive for grains or at least their icy 
mantles. This affects the statistics of the distributions and 
the various expectation values evaluated in this paper. For 
instance the velocity dispersion (Eq. (8)) is independent 
of the external scale and the density weighted velocity 
dispersion (13) is as well for any value (3 > 3, whereas 
it generally will not be for a fractal molecular cloud. In 
this case one should properly calculate the characteristic 
time it takes for grains to circulate into exposed regions 
from dense cold parts of the cloud core where the grain 
is shielded. The time of survival should then depend on 
the fractal dimensions of the hospitable and inhospitable 
parts of the cloud. One could think of this as consider- 
ably reducing the Reynolds number because in definition 
(6), the value for L should be the size of the high density 
(hospitable) regions. This means that for a fractal cloud L 
and the Reynolds number itself become stochastic quan- 
tities, for which expectation values must be calculated. If 
the regions with sufficiently high density occupy relatively 
small regions of space, e.g. if these regions are filamentary 
in structure, this reduces L on average and therefore also 
R e and thus T ioi . However the Reynolds number cannot 
be as low as to be of order unity on average over a molec- 
ular cloud because that would imply that on average the 
flow in a cloud complex is laminar which is clearly not the 
case. Therefore the estimate of Boland & de Jong (1982) of 
Ttot = L/a v for the residence time of grains in molecular 
clouds is too low. 
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3. Conclusions 

It is shown in this paper that t = L/a v severely under- 
estimates the actual time spent by grains in turbulent 
molecular clouds. In the limit of large Reynolds number 
for the turbulence, in clouds in which the velocity disper- 
sion deduced from line widths represents well the velocity 
distribution throughout the cloud, the residence time is 
Ttot ~ Q.35Re^ 4 L/a v . Unfortunately the Reynolds number 
is not a directly observable quantity. Using the definition 
(6) one can make a crude estimate however. The dissi- 
pation scale Id is determined by the microphysics of the 
mechanism responsible for the dissipation of the turbu- 
lence such as viscosity or magnetic diffusivity. Such mech- 
anisms usually become effective on length scales compa- 
rable to the mean free path of gas atoms or molecules. On 
larger scales the turbulence is essentially free and nearly 
dissipationless. Even at the low densities and low ioniza- 
tion fractions of cold interstellar matter this mean free 
path is not likely to be larger than 0.1 pc, at which scale 
ambipolar diffusion starts affecting the propagation of hy- 
dromagnetic waves (Mouschovias, 1991). The size L of a 
typical molecular cloud is more than 10 pc. Realistic val- 
ues of the Reynolds number for molecular clouds can thus 
be in excess of 10 3 , and have even been estimated to be 
of the order of 10 6 (Scalo, 1987), which implies that the 
residence time of grains can be more than 10 2 or even 10 4 
times longer than previously estimated. 

Taking as an example the clumps in the Rosetta molec- 
ular cloud which have a velocity dispersion of typically 
~ 1.5 km /s and a size of typically <~ 1.5 pc (cf. Williams 
et al., 1995) the time scale, for a Reynolds number of 10 4 , 
becomes ~ 3 10 8 yr. This is already of the order of the life 
time of a cloud, so grains formed in the centre of clouds 
like this will barely manage to migrate out of the cloud, 
let alone out of clouds that have a larger Reynolds num- 
ber. One might conclude that therefore molecules caught 
or formed on grains will not be returned to the ISM at all, 
since grains will not reach regions where photodesorption 
can play a role. However other effects might intervene. 
Vigorous turbulence is likely to lead to a higher frequency 
of grain-grain collisions, and also leads to heating because 
of dissipation of the turbulent energy. Inhomogcncity of 
a cloud produces a higher UV radiation field within the 
cloud than the cloud would have if it were homogeneous. 
What should be clear from this paper is that all physical 



and chemical processes in which grains are involved may 
well be modified by the presence of turbulence, and most 
certainly will have a much longer time to operate than 
estimated previously. 
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